Dynamical theory of superfluidity in one dimension 
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A theory accounting for the dynamical aspects of the superfluid response of one dimensional (ID) 
quantum fluids is reported. In long ID systems the onset of superfluidity is related to the dynamical 
suppression of quantum phase slips at low temperatures. The effect of this suppression as a function 
of frequency and temperature is discussed within the framework of the relevant correlation function 
that is accessible experimentally, namely the momentum response function. Application of these 
results to the understanding of the superfluid properties of helium confined in nanometer-size pores, 
edge dislocations in solid 4 He, and ultra-cold atomic gases is also briefly discussed. 
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Superfluidity and superconductivity are often associ- 
ated with the existence of long range phase coherence 
in quantum fluids. Nevertheless, long range phase co- 
herence is not a necessary condition for superflow, as 
the observation of superfluid response in two dimensions 
(2D) [1] demonstrates. In 2D 4 Hc films, torsional oscil- 
lator (TO) experiments have established [l[ existence of 
superfluidity (observed as a change in the resonance fre- 
quency of the oscillator), despite lack of the long-range 
order. The phase transition to superfluid phase is related 
to the binding of (free ranging) vortices and anti- vortices 
into pairs.as described by Berezinskii, Kosterlitz, and 
Thouless [2j. The 2D superfluidity without long-range 
off-diagonal order can still be understood in terms of 
the helicity modulus Q, which is a thermodynamic (i.e. 
static) property. However, superfluidity manifests itself 
experimentally as a dynamical property and in 2D, dy- 
namical corrections to the helicity modulus are im- 
portant in understanding experimental observations. In 
one dimension (ID), the dynamical aspect is even more 
important, since the helicity modulus vanishes altogether 
in the thermodynamic limit 31|. That is, dynamical ef- 



fects are not just corrections to the static picture, but 
are key to the understanding of superfluidity in ID. 

Recent TO experiments have detected superfluidity in 
long (0.2 — 0.5 /im) nanometer-sized pores filled with liq- 
uid He [5j, |6[ , where a suppression of the superfluid on- 
set temperature by pressurization and reduction of the 
pore diameter was observed. In optical lattices, it was 
found that a Bose-Einstein condensate of ultracold 87 Rb 
atoms exhibits coherent current oscillations 0, IH ■ How- 
ever, when confined to ID, the motion of the same ultra- 
cold degenerate gas becomes strongly damped even in the 
presence of a relatively weak periodic potential [|| . In the 
case of supersolid 4 He [§| it has been suggested theoreti- 
cally [HI U| and experimentally 12 1 that one likely ex- 
planation for the observations is related to the superfluid 
properties of edge dislocations in solid helium, which, as 



shown by Quantum Monte Carlo simulation [10j , behave 
as ID quantum fluids [l3l4l5j . These observations calls 
for a careful analysis of the notion of superfluidity in ID, 
despite the absence of helicity modulus in the thermo- 
dynamic limit. The superfluidity in ID is an essentially 
dynamical phenomenon, which also reflects peculiarities 
of dynamics in ID. Indeed, compared to 2D and 3D, dy- 
namics in ID tends to be much more constrained by the 
existence of conserved quantities. Recently, this has been 
shown to prevent complete thermalization [l6| or the to- 
tal decay of a current [17] in ID integrable systems. 

In higher dimensions, decay of superflow is caused 
by the motion of quantized vortices perpendicular to 
the direction of the flow. In ID, such a phenomenon 
corresponds to the creation of a topological excitation, 
namely a phase slip (PS), which 'unwinds' the phase 
difference imposed upon the system and whose impor- 
tance for the dynamical aspects of superfluidity in ID 
has been pointed out by several authors. In the frame- 
work of Landau-Ginzburg theory, a first calculation of 
the thermal production rate of PS was given in Ref. [l8[ . 
Later, these calculations have been extended to the quan- 
tum regime 11, [l9j]. In homogenous systems, the PS 
production rate is exponentially small at low tempera- 
tures, implying that the lifetime of the superflow in ID 
can be astronomically long. However, understanding of 
the suppression of superfluidity in the experiments men- 
tioned above 0, Q would require a finite PS production 
rate even at low temperatures. Moreover, the connection 
of the PS production rate to the experimental signatures 
of superfluidity, such as the response of a TO, remains 
obscure. 

In this Letter, we develop a theory of superfluidity in 
ID which emphasizes the dynamical aspects, as exper- 
imentally superfluid properties are probed at finite fre- 
quencies and are related to the dynamical momentum 
response function. To compute the momentum response 
we used the memory matrix formalism, which allows for 
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a perturbative treatment in the operators describing the 
quantum PS. By analyzing the momentum response as- 
suming a periodic potential (which is a relevant model 
for the 4 He systems of Refs. [1, [To! and the ID ultra- 
cold atomic gases in optical lattice|8l[20j). we show that 
the superfluid onset temperature decreases with decreas- 
ing the probe frequency (cf. Fig. [J) or decreasing the 
compressibility of the fluid (cf. Fig. [2J. The latter can 
provide an explanation for the pressure-dependent sup- 
pression of superfluidity observed in the nanopore exper- 
iments of Ref. [6| . 

Let us recall the description of superfluidity in terms 
of the momentum response function [2lj ]. which is the 
response of the system to the motion of the container 
walls. In higher dimensions, the fluid- wall interaction 
affects only the atoms in the neighborhood of the con- 
tainer, and it can be replaced by an appropriate bound- 
ary condition on the wall velocity field [2l| . The fraction 
of the fluid that is dragged along by a slowly moving 
container is given by the transverse part of the static 
momentum response function (2ll | which is defined as 
follows: Let II(r) = £ (r)W(r) - V^(r)V(r)] 
be the momentum current operator and X^v{ r T^) = 
-i»- 1 *(t)(|ll /i (r,t),n w (0,0)]> (ji,v = x,y,z) its re- 
sponse function, whose Fourier transform is a rank- 
2 tensor and for an isotropic fluid w) = 

(<W ~ ^)XT(q,u}) + ^XL(q,u), where Xt(l)(<7,^) 
is the transverse (longitudinal) momentum current re- 
sponse. The normal component density is p n = 
- 1 g-lim 9 ^olim w ^oXT(<7,^) (HI, Hi] (M is the particle 
mass). In contrast, in ID, the container wall affects the 
entire fluid. Therefore, its effect cannot be replaced by 
a boundary condition and has to be explicitly accounted 
for in the calculation of the momentum response. In fact, 
in ID, x(q> w ) i s a scalar which prevents the separation 
in a transverse and a longitudinal part. 

With an eye on the experiments 0, 



lflllj, we shall 



take a periodic potential to represent the container wall. 
In the experiments of Ref. [6] , the walls of the pore are 
covered by an inert layer of solid helium, which may 
be regarded as periodic (allowing for a disorder poten- 
tial is straightforward [15| and will not alter our con- 
clusions substantially). Furthermore, since we are in- 
terested in the low-temperature transport properties, we 
shall rely upon the Tomonaga-Luttinger liquid (TLL) de- 
scription of ID fluids [13" 15, 23, 24 1 where the low tem- 
perature/frequency degrees of freedom of the system are 
described by two collective (canonically conjugate) fields, 
9(x, t) and d x (j)(x, t) / 7r, which account for phase and den- 
sity fluctuations, respectively. The effective Hamiltonian 
takes the form H = Hq + H lrr , where 



#o = 7T I dx 
Ztt 



K{d x 0{x)) 2 + K- 1 (d x 4>(x)Y 



(1) 
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FIG. 1: (color online) Superfluid response as the frequency 
(wo = 2kHz 0). We used two terms for H QPS (cf. Eq. ©) 
setting Afcio = COOTa^ 1 and Afcu/Afcoi = 0.7. xo = 
M 2 vK/irh. The inset shows the imaginary parts of x( LJ ) 
(lower part) and the u-dependence of the dissipation peak 
temperature on a log-log scale. For these parameters, the 
dependence is (roughly) a power-law of ui following from the 
power-law prefactor of the memory matrix (see supplemen- 
tary material). Linear density and cut-off were chosen so as 
to conform to the experimental situation in 6], and we have 
taken the velocity of sound v = 200m/s and the Luttinger 
parameter K =8.1, resulting in onsets of superfluidity that 
agree well with experimental observations. 



= locity v. Using the TLL Hamiltonian ([TJ, wc find that 
x(w) = lim g _j. x(lt w ) — 0, implying there is no normal 
component and the system behaves as a perfect super- 
fluid at any frequency and temperature. Physically, this 
is because the Hamiltonian of Eq. (fTJ) neglects the exis- 
tence of quantum PS. This is, of course, a special prop- 
erty of (fT]), and the actual Hamiltonian of a ID fluid 
involves an infinite number of irrelevant (in the renor- 
malization group sense) operators H lrr and for the dis- 
cussion that follows we shall specialize to a particular 
subset of them, yielding the leading corrections to Hq'- 



H, 



QPS 



n>0,m 



Hvg n 



da; cos (2n(f>(x) + 2Ak nm x) 



Hq describes the properties of the system at T = 0, 
a ID fluid of compressibility K/(tvKVpo) and sound ve- 



. (2) 

which describes the effect of quantum PS, responsible 
for the decay of the momentum current; g mn are dimen- 
sionlcss couplings related to the strength of the periodic 
potential and the interatomic interactions; a ~ 1 is 
a short-distance cut-off; hAk mn = (2hnpQ — 2mG)h are 
the set of all possible (lattice) momenta carried by the 
PS (po being the fluid's linear density). As mentioned 
above, ID fluid is assumed to move in a periodic back- 
ground characterized by a minimum wave number G and 
the smallest | Ak mn | provides us with a measure of the in- 
commensurability between the ID fluid density and the 
wall potential 13J, [l5(- For Galilean invariant systems, 
G = and g n ,m^o = and vK — vf — hirpo/M. Irrele- 
vant terms like a^ n+m ^ J dx (d x <j)) 2n (d x 6) 2m (m+n > 2), 
etc., accounting for the curvature of the phonon disper- 
sion, do not, to leading order, contribute to the decay of 
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the momentum current. However, the leading irrelevant 
correction, H' irr = J dx (d x <f>) (d x 8) 2 Q, will be 

taken into account below when obtaining the low-energy 
form of momentum operator, II. 

In order to obtain the expression of the momen- 
tum operator at low temperatures/frequencies, we make 
a (time-dependent) unitary transformation to a frame 
where the walls are at rest 20]. This renders the cal- 
culation of the momentum response akin to the re- 
sponse of the system to an external gauge field pro- 
portional to the velocity of the walls v(t) [2(| [2lj j . 
Thus, U(x,t) — Mj(x,t), where j(x,t) is the parti- 
cle current operator. From the continuity equation, 
d±p(x,t) + d x j(x,t), and p(x) ~ po + ±d x <p(x) [H- 
1 1 51 ] - to leading order, j(x,t) = — ~d t (f>(x,t) = 
(«7r/i) _1 [Ho + Hi rr , <j)(x, t)] ]15|. Hence, including the 
leading order correction, n ~ — — J dx dt4>(x,t) — 



Jdx 
J = 



MvK 



AlvK 



d x <p(x) 



= J + —P, where 

VF 1 



j d x 6{x) 

— f dxd x 0(x) is the total particle (mass) cur- 
rent and P = ^ f dx d x <f>(x)d y 0(x) total energy current. 
Note that II is a linear combination of J and P, and 
these two operators are independently conserved by Hq 
(i.e. [H , J] = [Hq, P]=0 and (JP) = 0, where (. . .) 
is taken with respect to Hq). However, since both J and 
P do not commute with Hqps, in the presence of PS 
they become dynamically coupled and will acquire differ- 
ent decay rates. These effects can be taken into account 
within the memory matrix formalism [26l - l28j , which has 
been successfully employed to compute the AC conduc- 
tivity a (to) of charged ID systems [27L lijj]. 

In terms of the memory matrix M(u;T) the momen- 
tum response can be written as 

T) = Tr [V [ul + iM(oj; T)}' 1 iM{u; T)x(T)} 

(3) 

where V tj = (vK/v F y +: >- 2 = 1,2) and *(T) is the 
matrix of static susceptibilities (see supplementary mate- 
rial for detailed definitions). In Fig.[T]we have plotted the 
real and imaginary parts of the momentum response x( w ) 
against the absolute temperature, for different values of 
the probe frequency. In the inset we show the dissipation 
peak positions as a function of the probe frequency. The 
parameters of the system (see caption for more details) 
are chosen so as to reproduce onset temperatures com- 
parable to those experimentally observed in liquid He 
filled nanopores of Ref. @ when the probe frequency 
equals 2kHz [6]. As the probe frequency is decreased 
(corresponding to darker colored curves) , the onset tem- 
perature decreases. Indeed, this behavior can be antici- 
pated by taking the limit of ui — > + in (j3|), which yields 

X (w -> 0,T) = Tr [V X (T)\ = - (^)' 



Tr(k B T 2 ) 
Ghv 3 '■ 



which is in stark contrast with the vanishing result ob- 
tained by neglecting the PS. The limiting behavior at 
uj — > is also consistent with the vanishing helicity mod- 



ulus, namely absence of superfluidity in static sense. On 
the other hand, highly constrained dynamics in ID leads 
to superfluidity observable even at very low frequency 
such as 2kHz in Ref. 

In Fig. [5] we show real and imaginary parts of the mo- 
mentum response for several values of the TLL parame- 
ter K, which determines the compressibility of the fluid. 
The onset temperature is suppressed as the compress- 
ibility decreases (i.e. as strength of the the atom-atom 
increases). This is in agreement with the expectation 
that strong interactions tend to suppress the superfluid 
response. In the experiment of Ref. [6( the value of K is 
expected to decrease as pressure is applied to the system. 
Thus, the results displayed in Fig. [5] are consistent with 
the experimental observation that the onset temperature 
is suppressed by pressurizing the sample. Note also that 
since two separate currents are taken into account in ([3]) , 
depending on the parameters used, there exists the pos- 
sibility of two dissipation peaks with comparable weight, 
as displayed in the insets of Fig. [2] 
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FIG. 2: (color online) Superfluidity response as a function 
of the the Luttinger parameter K, which is determines the 
compressibility of the ID fluid. K ranges from 3.2 (dark) to 
10.2(bright) (Afcio = O-OOIoq 1 , and Aku/Ak w = 0.7). The 
insets (K = 6.2) show the effect of a larger value of the PS 
momenta (Afcio = O^ag 1 and Afcn/Afcio = 0.7), resulting 
in the appearance of two dissipation peaks. 



In ultracold atomic systems, the momentum response 
can be probed as described in Ref. [2(j. Thus, here we 
restrict ourselves to discussing how the momentum re- 
sponse is probed by a torsional oscillator (TO). Our start- 
ing point is Newton's equation of motion for a TO cell 
filled with liquid 4 Hc: 



dt 



[L z (t) + C z {t)} = -Ktp(t) - #W + r ext (t), (4) 



where L z (t) = Iotp(t) is the angular momentum and Iq 
the moment of inertia of the empy TO, ip(t) is the ro- 
tation angle, ip(t) = d<p(t)/dt, k is the restoring torque 
per unit angle, rj the friction coefficient, and r ex t(i) the 
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external torque driving the TO. C z (t) is the angular mo- 
mentum of the normal component of the helium sample 
(which is dragged along with the TO). Quantum mechan- 
ically, C z (t) = J dr(z x r) • (II (r, t)). For low rotation 
frequencies, this quantity can be computed within linear 
response theory: (II^(r,t)> = - J2 U J dtdr' x^{r,r' ,t - 
t') (<p(t')z x r') where x M „(r, r', t; T) is the momentum re- 
sponse of the liquid 4 He in the TO cell at an absolute tem- 
perature T 0. Hence, C z (t) = - J dt Xn (t -t';T) ip(t') 
where X„(t;T) = J drdr' g^(r, r')x^( r > r '> *5 T ) 

(jt,v = x,y,z), and g^(r,r') = (z x r) p (z x r')^.The 
TO response is defined by the relationship ip{t) — 
- J dtxTo{t ~ t')T cxt (t'). Hence, from g]), 



Xto( w ) = ^ 2 [ J o - Xn(w; T)} + irjuj - k. 



(5) 



In Ref. [25J, the TO response has recently been dis- 
cussed on phenomenological grounds. The expression re- 
ported there is identical to our Eq.([5]) provided we iden- 
tify g(uj]T) = — cj 2 Xn(w;T) where q(io \ T) is the back 
action function introduced in Ref. [25|. From ([5]) we 
see that the moment of inertia of the empty TO, Iq is 
corrected by 5I n = — Rex„(w;T) and the friction coef- 
ficient, is corrected by Srj n = — wlm x„(aj;T). To re- 
late x„(u;;T) to the momentum response of a ID fluid 
computed above, we imagine a straight ID channel lo- 
cated at r, filled with liquid 4 He, and oriented along 
the unit vector d. The momentum flow of the normal 
component is (II(r, t)) oc d and the velocity of walls 
of the channel is d • (z x r) ip(t). For a typical sam- 
ple size and TO oscillator driving frequencies ( ~ 10 3 
Hz [6|) this velocity field varies very slowly on the scale 
of the length of the channel (~ 0.3 [im 0). Furthermore, 
we assume that finite-size effects can be neglected (32| . 
Therefore, the momentum response of the channel is 
given by the q = q • d — > limit of x^i^T) com- 
puted for an infinite ID system. However, besides the 
liquid 4 He filling the ID channel, there may be an ad- 
ditional contribution to_y n (w;T) from other sources (in 
the experiment of Ref. [6] these would correspond to liq- 
uid helium filling the cavities between the nanoporous 
pellets). We shall assume that, around the onset tem- 
perature, these contributions only provide a weakly tem- 
perature and frequency dependent background signal. 
Thus, Xn(^;T) — Gox(w;T) + const, where Go is a ge- 
ometrical factor that measures the relative weight of the 
ID channel network to the total response of the sample 
in the TO cell. Hence, the change in frequency of the 
TO is 6u{T) = cu(T) - w +oc ~ G [1 + Re X (u ; T)/ X o], 
where ujq = y/ k/Iq and is the empty TO frequency 
(we neglect the difference between the TO resonance fre- 
quency and loq and uj +oc = ^n/(I + G xo), where xo = 
— lim:r_ i .o lim w _>.o x( w ; T) — M 2 vK/irh. The change in 
the quality factor is 5Q{T) ~ —Go Imx(w;T). Note that 
Go accounts for the distribution of orientations of the 
channels within the TO cell. Indeed, within the TO cell 
some ID channels will be oriented perpendicular to flow 



of the walls (that is, d • (z x r) ~ 0) and will not con- 
tribute to the change in the moment of inertia. 

Finally, let us briefly discuss the relevance of our results 
for the observed supersolid behavior in solid 4 He 0, [l2| ■ 
It has been suggested that in samples consisting of sin- 
gle crystals an explanation of of the observed supersolid- 
ity 0, [lH are edge dislocations. Indeed, associated with 
the onset of superfluidity in this system, there is a promi- 
nent dissipation peak and a stiffening of the crystals. The 
latter is believed to be related of the pinning of disloca- 
tions by 3 He impurities [12]. At sufficiently low temper- 
atures, the pinned dislocations become straight and be- 
have as ID superfluid channels that can be described 
as Tomonaga-Luttinger liquids with K ~ 5 [33j |. The re- 
sults reported here are fully applicable to such systems 
and improve on earlier theoretical treatments [lOj, 
Experimentally, it may be also interesting to further in- 
vestigate the similarities in the TO response between the 
nanopore systems 0, Q and the supersolid behavior of 
4 He single crystals. 

To summarize, we have reported a dynamical theory 
of superfluidity in ID quantum fluids, showing that in 
ID superfluidity is essentially a dynamical phenomenon 
related to the suppression of quantum PS at low temper- 
atures. Our calculations go beyond previous theoretical 
treatments by computing the experimentally accessible 
dynamical momentum response of the ID fluid which 
has been obtained in a particular case using the mem- 
ory matrix formalism, allowing us to take into account 
the dynamical coupling between the particle and energy 
currents. We have also demonstrated the explicit rela- 
tion between this response function to the measurable 
parameters of the torsional oscillator. 
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Supplementary Material 

The helicity modulus 

The helicity modulus Q is defined as 
T(T) = linu^ +00 T(T,L) where Y(T, L) = 
L(9 2 F(^i,T)/V) \ v=0 , where F(<p,T,L) is the 
free energy computed with twisted boundary conditions 
on the field operator: W(x + L) = e i<p 'J'(x). For a ID 
system described by the TLL Hamiltonian Ho it holds 
that 



T(T,L) = T 1 



(0) 



■v + OO 



e ^(0,e- 2e °/ fcBr ) 

= 2mz is a Jacobi theta 



where 3 O,«) = E„=-oo Q" e" 
function, and $3(2:, q) = d ^2 ; T = hvK/n and 
to = hirvK/L, where vif is the phase stiffness at T = 0. 
T(L,T) displays an onset at T ~ e 
infinite system size limit 



L 1 . Hence in the 



T(T)= lim T(£,T) = 



(7) 



excluding the possibility of a non- vanishing static super- 
fluid density. For a pore of R — 2.8 nm and lenght 
L = 0.3 /um , we find that when we subtract the inert 
layer on the inside of the pores that does not contribute 
to the supcrflow [5], at most T onsct ~ 0.18 K, which seems 
too low when compared to the experimentally measured 
T onsct ~ 0.65 K at a pressure of about 0.9MPa [6j. 
The memory matrix 

As defined in the main text, the momentum response 
can be written in the memory matrix formalism [26T - [28j ] 
in the following way: 

X(w; T) = Tr [V [wl + »M(«; T)]' 1 »M(w; T)x(T)} , 

(8) 



where x(^) 



diag{x,/,7, Xpp(^)} 



diagj— M f^ K , — ^g/fJP } is the matrix of static suscep- 
tibilities and the memory matrix, M(u);T) computed to 
lowest order in Hi rr takes the form 



M(w;T) 



h 2 n 2 

4ir 2 a 



n>0,m 



D nm {u-T)U mnX - l {T) (9) 



where 



^ 2nMvK y 
2nMvK A 7. 



2nMvK A jL 



(10) 



and 



D nm (a;;T) = [F ron (u;;T)-F mn (u; = 0;T)]/(ju;) (11) 
where 



+00 



77 m „(o;) = -i/T 1 y y dte iut-iAk m „x Cn (x,t) 
2alsmUn 2 K) /2-/, /; 2 



LM ^_.^ + ^Afc„ m ) i _ n2 ^ 

Zj| ^_.^ :! ;Afc„ ro ) ;i _ n2x 



4nk B T 



(12) 



being O n (x,t) = e- 2 *"^*), an d C n (x,t) = 
([O* (x, t), O„(0, 0)) , and 2?(x, y) is the Euler Beta func- 
tion. 



